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Orthogonal - unitary and symplectic - unitary crossover ensembles of random matrices are relevant 
in many contexts, especially in the study of time reversal symmetry breaking in quantum chaotic 
systems. Using skew-orthogonal polynomials we show that the same generic form of n-level correla- 
tion functions are obtained for the Jacobi family of crossover ensembles, including the Laguerre and 
Gaussian cases. For large matrices we find universal forms of unfolded correlation functions when 
expressed in terms of a rescaled transition parameter with arbitrary initial level density. 
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Random matrices have found applications in a wide 
class of systems exhibiting broad range of behavior [H-O ■ 
The Gaussian ensembles of random Hermitian matrices 
have turned out to be particularly useful and have been 
studied extensively and applied successfully to many sys- 
tems. In recent years, however, non- Gaussian random 
matrix ensembles have drawn considerable amount of at- 
tention. The reason behind this is that, apart from the 
mathematical interest of extending the results known for 
Gaussian ensembles to the non-Gaussian ones, they turn 
out to be of great importance in studying physical sys- 
tems where these appear naturally and hence are more 
appropriate in modeling the system. For example the 
Laguerre (or Wishart) ensembles comprise non-negative 
definite matrices of the type A'' A (or ^yl^) where A is 
a rectangular Gaussian random matrix [HI, M- Commu- 
nication theory 0, amorphous systems [3, [l0| . and 
multivariate analysis of chaotic time series are a few 
examples where Laguerre ensembles of random matrices 
appear. Jacobi ensembles are generalizations of Gaussian 
and Laguerre ensembles and serve, inter alia, as models 
for transmission matrices in quantum transport problems 
0i [13 ■ Depending on the invariance under orthogonal, 
symplectic and unitary transformations the ensembles 
are referred to as orthogonal (OE), symplectic (SE) or 
unitary (UE), the first two being applicable to time rever- 
sal invariant (TRI) systems and the last to time reversal 
non-invariant systems [H-Il]- The corresponding Gaus- 
sian ensembles are abbreviated as GOE, GSE and GUE, 
the Laguerre ensembles as LOE, LSE and LUE, and the 
Jacobi ensembles as JOE, JSE and JUE. Our purpose in 
this work is to study non-Gaussian ensembles of random 
Hermitian matrices which interpolate between ensembles 
with orthogonal invariance and unitary invariance. The 
crossover is governed by a parameter r with r = for 
the orthogonal case and t = cx) for the unitary case. We 
similarly study symplectic - unitary crossover, again as 
a function of the parameter r. These ensembles serve 
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as models for systems with partial TRI breaking. In ap- 
plications to quantum chaotic systems r is a measure of 
TRI breaking, r = being the TRI case. There are also 
applications in other systems mentioned above j4j, .8-10]. 
The Gaussian ensembles with the same (OE-UE and SE- 
UE) crossovers have been studied earlier [l^ and applied 
to nucle a.r sp ectra [l^ and to spectra of quantum chaotic 
systems jl5| . See also [l^ for the same transitions in the 
circular ensembles. 

We consider here the problem of (OE-UE and SE-UE 
) crossovers in Laguerre and Jacobi ensembles. We find 
that the expression for the joint-probability density (jpd) 
of the eigenvalues for these transitions has a generic form 
applicable to all crossover ensembles belonging to the Ja- 
cobi family. We give exact compact expressions for the 
jpd of eigenvalues and n-level correlation functions for 
finite-dimensional matrices. This includes results for the 
skew-orthogonal polynomials appropriate to the transi- 
tions which are generalizations of those studied earlier 
[13] ■ We also propose a method for obtaining directly 
the large dimensionality limits of the unfolded two-point 
kernels and n-level correlation functions. The unfolded 
results are found to be the same as those in Gaussian 
and circular cases with the rescaled transition parameter 
A now appropriate to the Laguerre and Jacobi ensembles. 
We stress that the results are independent of the initial 
level density, as the latter affects only the unfolding func- 
tion of the spectrum and the rescaling of the parameter. 
Finally we discuss briefly applications of these ensembles 
to some of the systems mentioned above. 

We consider TV-dimensional Hermitian matrices with 
eigenvalues Xj where j = 1,...,N. The eigenvalues sat- 
isfy oo > Xj > — CX), Xj > 0, and 1 > Xj > —1 for 
Gaussian, Laguerre and Jacobi ensembles respectively. 
We define the r-dependent jpd, V{t) = V( x^ , x at; t), 
using Dyson's Brownian motion model @, [l8l - [20j . We 
have dV{T)/dT = —£V{t), where £ is the Fokker-Planck 
operator with equilibrium jpd Vcq = Vcq{xi, ...,xn)- 

— 1/2 

The similarity transformation ^ — "Poq P and "H = 

1/2 1/2 

■peq CVcq leads to the evolution equation of ^ involv- 
ing Calogero-Sutherland type Hamiltonians [21, 22] for 
the ensembles mentioned above 0] . The ground state of 
T-L is non-degenerate with zero energy so that Vcq exists. 
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For transitions to UE, the interaction term disappears 
and then H can be expressed as sum of N single-particle 
Hamiltonians Hx representing noninter acting fermions. 
Thus H = X^jLi '^^j —^0- Here £o is the ground state en- 
ergy of the iV-fermion system. The T-Lx are non-negative 
definite operators and, for the Gaussian, Laguerre and 
Jacobi ensembles, are respectively d, 0, [l^, l20| 



H'l — 
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where a>— 1,6>— lin® and ([3]). For example, 
2a + 1 = A^' — > in ([2]) for transitions correspond- 
ing to the Laguerre ensembles A'^ A, where ^ is iV' x 
dimensional. Here the matrix elements of ^ = A{t) are 
independent complex Gaussian variables with zero mean 
and same variance. For OE-UE transition the real and 
imaginary parts are also independent and have the vari- 
ances (l+e"^)/4and (l-e~'^)/4 respectively. For SE-UE 
transition, same properties are valid with quaternion real 
and quaternion imaginary parts. 

Using the Fokker-Planck operator £, we have, for any 
arbitrary initial jpd 'P(O), 

r{T) = e-^^V{0) = <^e-«-7'-i/2p(o). (4) 

This expression governs the transition from any arbitrary 
initial jpd to the jpd corresponding to UE as the final or 
equilibrium density and can be efficiently utilized to solve 
the problems of OE-UE and SE-UE transitions. In (g]) 
■peq corresponds to the jpd of UE, 
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Vcq OC IAatP Y\_ Wcq{Xj), 



(5) 



where Wcq is the square of the ground state wave function 
of Hx and Ajy — Ylj^^.{xj — Xk) is the Vandermonde 
determinant. For the three ensembles, Woq is given by 



Wcq{x) 



for GUE, 
for LUE, (6) 
(1 - a;)2«+i (1 + x)^^+^ for JUE. 



In the SE-UE transition N is even because the eigen- 
values are doubly-degenerate for t = (Kramers degen- 
eracy) . We consider here N as even for OE-UE transition 
also. The initial jpd of eigenvalues is 



N 

P(0) cx An PW°>{ )] J|u;(a;,;) 



(7) 



Here j. A; = 1, 2, TV, w{x) is the initial weight function 
with boundary conditions similar to those for Weg(x), 



^W(x,y) 



e{x — y) for OF, 
-5'{x - y) for SE, 



(8) 



2e(x) is the sign of x and Pf(i3) is the Pfaffian of antisym- 
metric matrix B. Then Pfaffian of C?'^"^ (x, y) along with 
Ajv leads to the factor of \An\ for OE, whereas for SE 
it gives rise to lAr/y/-?)!"^ a long with the delta functions 
for the degeneracy [ll. JlSl. fla ]. Using (4), (5) and the 



expansion of Pfaffians [l| , we get the jpd for arbitrary r, 

JV 



V{t) (X e^o^Ajv Pf[g^^Hxj,Xk)]l[w{x,). 



(9) 



Here ^('^^ is defined in terms of one-body operator Ox, 
g^-\x,y)^OxOyg^''\x,y), (10) 



_ Vweqjx) W{x) 
^x — / N 



w{x) 



(11) 



We now introduce the skew-orthogonal polynomials 
q'^^\x), the weighted polynomials (jJ'p {x) — w{x)qj^\x), 

and the integrated functions ipj^\x)- With j,k = 
0, 1, 2, we have 



(l)^f\x)'>p[''\x) dx = Zjk, 



i^\^\x)^ / g'^'\x,y)<l,^;\y)dy, 



(12) 



(13) 



where Zjk = —Zkj equals 1 for A: = j -|- 1 with j even, — 1 
for fc = j — 1 for j odd, and for \j — fc| 7^ 1. [The inte- 
grand in (IT^ should be antisymmetrized in the symplec- 
tic case when w{x) diverges.] Note that the polynomials 
are skew-orthogonal because ^/^"^^ is antisymmetric. For 
T = 0, ^ and (HSl) give back the OE and SE definitions 
[li[l3. We can show from ([H]) and ^ that 



00 

E [4l\^)4l\iiy) - 4lU^)4l\y)] , (14) 



ti=0 



and then, using (fTU|) and IT^ . we get 
i;^'\x) = O.^f (x), <l>^;\x) = (Ot)-Vf (x). (15) 

These functions can be calculated using the spectral de- 
composition of Hx- 

With special choices of the initial weight function we 
get compact expressions for (f>'^^\x) and 'il'j^\x), similar 
to those for r = [l^- The choice w^x) = Wa{x) = 
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x^e ^ is appropriate for the LOE-LUE transition. We 
find, with ipo{^) given by (|f 51) . 

2a+l/2 

a;2M 



4;^^) = e^''^^^Wa{x)L'i;+'\2x), (16) 



Then, using Dyson's theorems P, HI], the n-level corre- 
lation function can be expressed as a quaternion deter- 
minant (Qdet), 

i?„(a;i, ...,a;„;r) ^ Qdet[a'^'^\xj,Xk)]j,k=i,..,n, (24) 



g _ ^iV-Y^(2a;) 

as,, 1 V ; 

■— 5 V'2/.-2W, (M>0), (17) 



4;Vi(^) = e''^^^—^wa{x)\e^{2fi + 1) 
a2p L 

xL^l^+,'\2x)-e-^2^, + 2a+l)Li';+,'\2x) , (18) 



Here Lf"+^)(a;) are the associated Laguerre polynomi- 
als with weight function W2a+iix) and normalization 
f^(2a+i) ^ ^^^^2 ^ Y{j + 2a + 2)/r(j + 1). The choices 

w{x) = e~^^/2 and (1 — a;)''(l + a;)'' give compact re- 
sults respectively for GOE-GUE and JOE-JUE transi- 
tions, which, for t = 0, coincide with the results in [l7| . 
For GSE-GUE, LSE-LUE and JSE-JUE transitions, the 
choices u;(x) = e'^^'/^^ x^+^e""" and {1 - x)''+\l + xf+^ 
give compact answers respectively. 

We now evaluate the n-level correlation function i?„, 
defined by 



2a+3/2 

1 

Q!2/i 



,+i(:r)4^;+^^(2x). (19) 



dx. 



n+l 



...dxNPir), (20) 



{N -ny. 

where n = 1, 2, N. For this, we introduce the kernels 

i-r) 
N 



S^j^\x,y), A^^\x,y) and B^^\x,y) along with their r- 



V-^' f-*' ^^AT 

evolutions. We have 



(N/2)-l 



E ^2;^(^)4;Vi(2/) 

:(Ot)-iO,5(,°)(x,y), 



(21) 



where, for both transitions, a^'^\x,y) is given by 
aM(x,y) 



'S^^\x,y) A^^\x,y) 
B^^\x,y) S^j^\y,x) 



(25) 



For large N, the level density Ri{x-t) = 5*^^(2;, x) 
undergoes a smooth transition from r = to r = 
oo as a function of r (or Nt in some cases). How- 
ever the unfolded correlation functions i?„(ri, r„) = 
Rnixi, ...,Xn)/Riixi) ■ ■ ■ Ri{xn), where Xj = x + 
rj/Ri{x), undergo the transition for much smaller r, 



given by (|^ below. Consider r = {x — y)Ri{x) where 
Ri{x) = i?i(a;;0) is the level density corresponding to 
the initial ensemble. For large N we have Tlx — Hy « 
and Hx + Hy « 2f{x)Rl{x)d^ / dr^ , where f{x) = 
— 1/2, — a;/2, — (1 — x^) for Gaussian, Laguerre and Ja- 
cobi cases respectively and is the factor appearing with 
/dx^ term in T-L^- Let 



= ^-Tf{x)Ri(x). 



(26) 



For large N, s'^j^^ / Ri{x) in (1^ becomes independent of 
T and has the limit 



S{r) 



sm Trr 



1 



cos kr dk 



(27) 



7rr TT Jo 

for both the transitions. On the other hand 
/{Ri{x)f and sj^' have the limits 



A{r-X) 



^(r;0). 



S(r;A) =6^^ ^ /^'- i?(r;0) 



(28) 



(29) 



We use OE and SE results [U, [11, [T^] for A = 0. Thus, 
we find 

1 

A{r-\) = / fcsinfer e^^'*-'' dk, (30) 



(JV/2)-l 

-\Ol)-\Ol)-'A%\x,yl 



(22) 



M=(JV/2) 



(23) 



B{r;X) = — 
n 



sin fcr 



-2A^fc^ 



dk, (31) 



for OE-UE transitions. Similarly we obtain 



A{r:X) = 

TT 



1 r sin kr ^^2^2 



dk. 



B(r; A) = / fc sin fcr g-^^'*^' dfc, 

Jtt 



(32) 



(33) 
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for SE-UE transitions. For both transitions 

fi„(ri,...,r„;A) = Qdet [cr(rj -rfc;A)] .^^^ „ , (34) 



a(r;A) 



' Sir) A{r;X) 
i?(r;A) 5(r) 



(35) 



These results are same as those for the Gaussian and cir- 
cular ensembles [ij, ■ Note that the two- level cluster 
function is given by Y2{r) = 1 — i?2 = 5*^ — AB where 
r = ri — r2. The number variance and other two-point 
fluctuation measures derive from Y2 and have been used 
in the study of TRI breaking in complex nuclei and 
quantum chaotic systems [l^- We believe that similar 
transition will be obtained in the vibrational spectra of 
amorphous clusters [1, • In some of these applications 
Laguerre ensembles are a priori more appropriate. 

In the study of quantum transport in (mesoscopic) 
chaotic cavities, the transmission eigenvalues Tj after the 
transformation Xj — 2Tj — 1 are described by the Jacobi 
ensembles [1, [13] ■ For JOE-JUE transition the weight 
function changes from (1 + x)^ to (1 -|- x)'^^^^ , where 
2b + 1 = |iVi - A^2| and N = min(7Vi, A^^z)- Here Ni 
and N2 are the number of incoming and outgoing chan- 
nels. Again the (f> and ip functions can be found explic- 
itly using the above operator method and all quantities 
of relevance for conductance fluctuations Q can be de- 
rived as a function of the TRI breaking parameter t for 



arbitrary 7Vi,iV"2. For example, the variance of conduc 
tance is given by (1 -t- e-^'-^^+^^^''){NfNi/{Ni + N2f 
for large A'^i and N2- Our results agree with those in 
for T — 0, 00. We also remark that it should be possi- 
ble to derive a long-range two-point correlation function 
analogous to those given in ,24] for the transitions also. 
In fact such expansion has already been given for tran- 
sitions in the Gaussian ensembles [2], jl^]. Finally, our 
LOE-LUE transition results are directly applicable in the 
calculation of the Shannon capacity of a^multiple input 
multiple output communication channel [8] for arbitrary 
number of transmitting and receiving antennas. 

To conclude, we have shown that the method of skew- 
orthogonal polynomials is valid for a large class of OE-UE 
and SE-UE transitions. The polynomials can be calcu- 
lated from those for the initial ensembles in terms of one- 
body operators. Our large-iV results prove the univer- 
sality of spectral correlations in the crossover ensembles 
of the Gaussian, Laguerre and Jacobi types with arbi- 
trary initial level density in terms of a single parameter A. 
(Generalizations to arbitrary initial level density can also 
be done for the circular ensembles.) The finite- iV results 
are useful in communication theory 0, 1^ and mesoscopic 
quantum transport problems [3| , and the large- results 
are useful in quantum chaos studies El] • Details of 
this work appears elsewhere [2^] . 
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